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Abstract 

We have calculated the KK — > KK scattering amplitude to next 
to leading order in Chiral Perturbation Theory. Then, making use of a 
unitarization procedure with one or several coupled channels (mr^KK 
in our case) we have calculated the tttt —f tttt, tttt — > KK and KK — > 
KK S and P waves in good agreement with the experiment up to 
1.2 GeV. 

The TTTT scattering lengths with isospin and spin (I, J) equal to (0,0), 
(1,1) and (2,0) are also calculated in agreement with experiment and 
former Chiral Perturbation Theory calculations. 

Finally we have employed these amplitudes, making use of an Omnes 
representation, to calculate the scalar and the vector pion form factors, 
obtaining a good agreement with the available experimental data. 
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1. Introduction 



The low energy effective theory of the strong interactions (QCD) is Chi- 
ral Perturbation Theory (ChPT). The chiral symmetry constraints are 
a powerful tool, enough to determine the low energy matrix elements in a 
systematic way ^, ^ |5|, ^. In this way one can evaluate the scatter- 
ing amplitudes corresponding to the lightest octet of pseudoscalar mesons 
(tt, K, rj) in a loop expansion. This has already been done for the tttt — > vrvr 
scattering amplitude for the SU{2) case in |^] and extended to SU{3) in |^. 
The Kit — > Kit scattering amplitude is also evaluated in Q. 

This loop expansion can also be seen as an expansion in powers of the 
masses and external tetramomenta of the pseudoscalars over some typical 
scale around 1 GeV, so that the series is, in principle, only useful for low 
energy. Furthermore, in the meson-meson scattering, it is common the ap- 
pearance of resonances with a typical mass around 1 GeV, which obviously 
can not be reproduced in a power expansion because they correspond to 
poles in the T matrix. The region of validity of the series depends on the 
channel (I, J) considered, where I accounts for the isospin and J for the an- 
gular momentum. A way to get a hint about this convergence problem is 
to calculate the next to leading contribution {0{p^)) and compare it with 
the lowest order contribution(C'(p^)). This seems logical and is the rea- 
son advocated in |^ to justify their Kit 0{p^) calculation even though the 
threshold for this reaction is 0.64 GeV and typically the chiral expansions 
break around ^/s ~ 0.5 — 0.6 GeV. 

In this paper we first calculate the KK KK scattering amplitude 
at 0{p'^) in ChPT. This implies to evaluate two amplitudes due to the fact 
that both 1=0,1 are possible. In this way we calculate the K+K' K+K' 
and K^K'' K^K^ scattering amplitudes. It is not necessary to calculate 
the K^K^ K^K^ amplitude because in the isospin limit T{K^K~ — > 
K+K-) = T{K^K^ ^ K^K^). 

This calculation seems in principle to be unreasonable since the KK 
threshold is almost 1 GeV, and furthermore, for 1=0,1 at this energy the 
/o(980), ao(980) resonances appear which couple strongly to the KK chan- 
nel. In fact, after making this Oij/^) calculation, we will see that, at this 
energy, the 0{p^) result is larger than the 0{p'^) lowest order calculation. 
Obviously, what all this is telling us is that one cannot rely in the per- 
turbative chiral expansion at these energies. This implies that one cannot 
compare the predictions directly with the experiment, using the 0{p^) plus 
©(Z) ChPT amplitudes. 
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In the last years, there have been several attempts to extend ChPT to 
higher energies. One of them uses an effective field theory with explicit 
resonance fields as degrees of freedom |, g. This resonance chiral effec- 
tive theory has one drawback: applying only the constraints coming from 
the symmetry the number of free parameters grows so fast with higher or- 
ders that in practice predictions become impossible. However, using this 
lagrangian at 0{p^) and imposing some short distance QCD constraints 1^] 
one can obtain interesting and good results. For example, the study of the 
vector pion form factor done in [^]. 

Another attempt (the one we are going to use here) is to unitarize the 
ChPT amplitudes. In the Inverse Amplitude Method (lAM) was 

developed and used to study the tttt and Ktt scattering. This method has 
the problem that only elastic unitarization is employed, so that, for those 
cases where inelastic thresholds are important, as the KK one in / = J = 0, 
the improvement over ChPT is poor. But for those other channels predomi- 
nantly elastic, as the I = J = 1, the improvement is remarkable and, in fact, 
in this channel the p resonance is obtained. In order to use the lAM one 
needs the 0{p^) ChPT amplitude. More recently, in [^4|, a resummation of 
the chiral series was done inspired in the Bethe-Salpeter (BS) equations for 
the J = sector making use only of the 0{p'^) ChPT amplitudes and im- 
posing unitarity with coupled channels, vrvr and KK. The agreement with 
the I = J = and / = 1, J = experimental data was very good, re- 
producing the presence of the /o(980) and ao(980) resonances respectively. 
However, with this method the higher orders are obtained through unitar- 
ity loops in the s channel which are subleading in the large Nc limit and 
then fails in the vector sector where the large Nc limit works very well. 
The leading part in this latter limit are tree level contributions appearing 
through the chiral lagrangians. In fact in |^ it is seen that in the J = 1 Ktt 
scattering the 0{p^) corrections are clearly dominated by the polynomial 
counterterms contribution coming from the 0{p'^) chiral lagrangian. This 
same phenomenon is also seen in the P-wave vrvr scattering and in the vector 
form factor ||2|, 15, 17|. It is clear then that a unitarization proce- 



dure which could include the success of the lAM and BS approaches, or 
equivalently, a method that could handle unitarization in coupled channels 
incorporating both the 0{p'^) and 0{p^) chiral amplitudes would be very 
welcome. This in fact has been recently done in |18| and the resonances 
appearing in both J = 1 {p{770) and K*{890)) and J = (ao(980), /o(980)) 
channels were reproduced. This is in fact the method we are going to use 
here. The novelties are that while in the 0{p'^) chiral amplitudes were 
approximated in a way inspired in 1 14 1 , here we are going to use the exact vrvr 
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and Ktt scattering amplitudes [Q, |^ and the KK ones, which we calculate 
here to 0{p^). With these ingredients we obtain the scattering amplitudes 
for the channels (I,J)=(0,0), (1,1) and (2,0), generating dynamically the /o 
and p resonances. The agreement with experiment is good up to around 
^/s ~ 1.2 GeV as it will show. 

Finally, making use of the calculated vrvr vrvr phase shifts we derive 
an Omnes representation M, EOl for the scalar and vector pion form factors 



based on the Watson final state theorem |21|. For the case of the vector 
form factor experimental data are available and, with our calculation, the 
agreement is rather good. In the case of the scalar form factor we compare 
the calculation making the unitarization with and without including the KK 
channel and the differences are very significative, even at energies around 
|i/s| = 500 MeV. The /o(980) resonance appears clearly in the scalar form 
factor of the pion. 

2. The KK —>■ KK scattering amplitude to next to leading order 
in ChPT 

As it was stated before, in order to calculate the KK amplitude, we need 
to evaluate two amplitudes. These two independent isospin amplitudes can- 
not be connected by crossing symmetry because they have different absolute 
values for the strangeness. This is contrary to what happens in Ktt scatter- 
ing with / = 3/2, / = 1/2. 

We calculate the amplitudes K+K' K+K' and K+K' K^K^ 
which we denote by Tec and Ten respectively. 

The scattering amplitudes with definite isospin T^^^ can be written in 
terms of Tec and Ten in the following way: 



r(°)(s,t,M) = Tce{s,t,u) +Ten{s,t,u) 

TW{s,t,u)= Tce{s,t,u) -Tcn{s,t,u) (1) 

We now proceed to describe the calculation scheme for these amplitudes up 
to 0(/). 

At lowest order one has the ChPT lagrangian at 0{p'^) 



c^ = !l(d^u^d^'u + m{u + u^)) (2) 



4 

where () stands for the trace of the 3x3 matrices built from C/($) and 7W, 
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[/($) = exp 



iV2 



where $ is expressed in terms of the Goldstone boson fields as 



7r~ 



The mass matrix M is given by 





-ml 





M = 





ml 




. 










2rn|^ — m: 



K+ 

2 



(3) 



(4) 



(5) 



in the isospin limit. Where m means bare masses. 

From this lagrangian we can evaluate the lowest order contribution to 
the KK scattering amplitudes 



Tcc,2is,t,u) 
Tcn,2is,t,u) 



2.2^4 2 

-m^ + -m^ - u 
-m^ + -mK - u 



(6) 



where the subindex 2 means 0{p'^). 

At fQ = f^ = 93.3 MeV and rfiK = mx = 495.7 MeV. But when 

we go to next order these equalities do not hold. This is the reason why we 
keep the distinction between bare and physical masses. 






rji I rjiVJ T 

Figure 1: Diagrams at 0{p^) 

At 0{p^) one has to calculate the diagrams schematically shown in Fig. 

1. TJ represents contributions coming from the £2 ChPT lagrangian with 
six fields and a tadpole loop. The represents the loops constructed from 
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the £2 amplitudes with four fields appearing in the vertices of the loop. We 
will call this contribution unitarity loops because it makes the amplitude 
unitary at 0{p^). These loops include contributions from loops in the s,t 
and u channels, as shown in Fig. 2. 




K K K K K K 

s-channel t-channel u-channel 

Figure 2: Diagrams for the s, t and u channels 

Finally the amounts for the 0{p'^) polynomial contribution coming 
for the £4 ChPT lagrangian, which can be written as, 

£4 = Li (^df^u^d^uy + L2 (d^u^d.u) ld>'u^d''u) (7) 

+L3 (d^U^d''UdyU^d''u) + L4 (a^C/ta^f/) (u^M + M^u) 
+L5 (df,U^d>'U {u^M + M^U^ ) + Le (u^M + M^U^f 

+L7 iu^M - M^Uy + Lg {m^UM^U + U^MU^m) 

When taking into account the wave function renormalization, Fig. 3, and 
the relation between bare and physical masses and decay constants, other 
0{p'^) contributions come from the lowest order amplitudes. 




Figure 3: Wave function renormalization 

The relation between mx and rhx and the one between /o and /tt at 
0{p^) can be obtained from 
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The final amplitude up to next to leading order is then obtained summing 
all the former contributions. We express it divided in four parts: T2, TJ', 

and T^" where the subindices indicate the order in powers of momentum. 

In TJ' and we have also included, in addition to the one coming from 
Fig. 1, the 0{p^) contributions from the renormalization of wave function, 
masses and decay constants. The ones with Lj parameters are included in 

and the rest, of tadpole type, in Tj^ . 

The amplitude for K~^K~ K^K'^ is: 

Tcc,2{s,t,u) = (8) 



a: 



288/47r2 



r(8m|: + m^-3u) + 



A 



K 



288/47r2 



(8m|: + 3u) 



Qu) 



(9) 



-20m% - 2s^ - 2su + u{3ml + u) + mj^{8s + llu - Ami) 



Al 



+ 



[ kis) 
1 1536/47r2 



8m|-(s - Ami) + - 4it) + 8m^(s + 2u) 



^^^^^ 5(8m^ + s{s -u)- Amj^t) 



384/47r2 
B'nis) 



+ 



4608/47r2 



{8m\ - 9s)' 



768/^2 



B^iu) 



32/>2 



{u - 2m\^Y (10) 



= ^1 74 + s' + + 4m|^n) 

J TT 

+L2^(s2 + ^2 + - Am\u) 
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4 
'74 " 

J TV 



L3^{-8m% + 8"^ + f + Amj^u) 
i\u 8m\u 

fl — Jl l-(2-L6 + -o/ .4 

/tt /tt /tt 



-La 



32m 



(11) 



The amplitude for K+R- K^K^ is 



2m|^ - n 

J-r.n..2 — 



(12) 



576/47r 



288/>2 



(lOmlf - 3t) 



+ 



A''' 

(20m|; - - 6u) 



576/47r2' 



(13) 



rpU 

-'■cn,4: 



24m% - 6mls - - 2su + + 2m'j^{iml - 7t) 



384/47r2 



a: 



1152/47r2 



:{38mj^ - 2ml - 6s - 15u) 



576/47r2 



(10m|f + 3s - 6u) + 



1152/47r2 



(22m|: + 2m^ - IS-u) 



(8m|f (s - 4m2) + s(7s - 4i) + 8ml{s + 2i)) 



1536/^2 

^^^{s -u){t- 4ml) + + sis -u) + 4mlis + u)) 

(-8m|, + t{t -u) + 4mlit + u)) + ^^{u - 2mlf 



384/47r2 



64/>- 



4608/^2 



(8<-9-)'-JiS(3--4ml,)^ 



768/47r- 



384/47r2' 



(14) 



cn,4 = -^174(5 -2m|f)2 

/tt 

+L2^(-8m^ + + '"^ + 4mlfs) 
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+L3 — {-8m% + s^ + t^ + Am\u) 

Jtt 

+L, — {-2Ami + Umls) - L5— ^ 



+(2L6 + L8) 



16m^ 



#4 



In the above formulas we have used the quantities 



(15) 



Ap 



-nip 



-1 + ln 



rriT 



(16) 



A^/^(g,m^,m^) / m| + - g + A^/^(g, mj,, m^) ' 
^"^^^^ 2s \ml + ml-s-\y^s,ml,ml) 



+2 -In 



m} 



+ 



In the equal mass limit ( p!7|) reduces to 

2 ^ 



y^rrip + uIq — 
ml-m'^Q + s / ml 
\ ml, 



B'p{s) = 2 - In ( ^ - a{ml, s) In 



a{m'p, s) + 1 
a{m'p, s) — 1 



(17) 



(18) 



where 



A(s, mp,mQ) = [s - {mp + rnqf] [s - {mp - mq)^] 



a{m'p,s) =dl--^ (19) 



The functions (16) and ( p^ ) come from the Passarino-Veltman integrals with 
one and two propagators [p^]. 

It is interesting to note that Lj does not appear in T^. This also happens 
in vrvr — > vrvr and Kit — > Kn 0{p'^) scattering amplitudes. Lq and Lg appear 
in the combination 2Lq + Lg as a consequence of the Kaplan and Manohar 



symmetry [23 



As it was stated before, the corrections coming from the 0{p'^) calcula- 
tion are, at least, as large as the lowest order contribution itself. For example 
in the KK threshold: 



56.5 
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Tcc,i = 73.6 + ii 36.74 



Unambiguously this means that a perturbative calculation is useless in this 
region and that some non perturbative scheme should be used in order to 
compare with the experimental phenomenology. 

3. Unitarization of the vrvr and KK amplitudes 



We already stated in the introduction that we are going to use a uni- 
tarization procedure recently developed in and thoroughly used and 
applied to phenomenology in |^^, where the S and P waves meson- meson 
amplitudes were reproduced successfully. However, in both works, the 0{p^) 
amplitudes were approximated. Here we are going to use our calculated 
0{p^) KK amplitude together with the vrvr and Kir 0{p^) scattering am- 
plitudes given in Q. 

In a rather general scheme was derived to obtain final unitarized 
amplitudes from the C(p^) and C(p^) ChPT scattering amplitudes. For our 
case with two channels, labelled as 1 for vrvr and 2 for KK, the amplitude 
with definite I, J is given by 



^2 ~ -'4 



rf'*^^ (20) 



where the different amplitudes are actually 2x2 matrices for 1=0, 1, ac- 
cording to the above labelling of the states, and just numbers for 1=2. 

Note that (^) was obtained from an expansion of the l/T matrix ampli- 
tude, which has a zero when T has a pole. In this way one expects that the 
l/T expansion in powers of masses and momenta is meaningful even above 
the resonance region, where the low energy expansion of T has no sense 

The projection in definite angular momentum is given by 

T^'''^ = ^ j\ d{cose)T('\s,t)Pj{cos9) (21) 

where = 1 for KK KK, iV = ^/2 for vrvr ^ KK and iV = 2 for 
vrvr vrvr, since in the isospin formalism the pions are identical particles. 
With this normalization, in our case unitarity reads, for (I,J) = (0,0) 



4ml < s < 4m\ : 
4:m\ < s < 4m^ : 



ImTii = a{ml,s) \Tii\ 



(22) 



Imri2 = cr{ml, s) TuT^^ + f^("^i^> ■?) ^12^22 
Imr22 = a{ml,s) \Ti2\'^ + a{mj^, s) \T22\'^ 
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For (I,J) = (1,1) 



ImTii = a{ml,s) \Tu\'^ + a{mK,s) \Ti2\'^ (23) 



Ami < S < '^ITT'K 

'im\ < s : 



And for (I,J)=(2,0) 



ImTii 
Imri2 
Imr22 
ImTii 



cr(m^,s) |Tii| 



(24) 



a{ml, s) TuT^2 + (^i^K^ «) T12T22 
a{ml, s) |Ti2|^ + cr{mj^, s) \T22\'^ 
aiml,s)\Tu\^ + a{ml,s)\Ti2\^ (25) 



4m^ < s 



ImTii = a{ml,s) \Tii\ 



with fj given by (|T9|). 

In matrix notation ( p3|) and (25) become 

ImT = TaT* 
with a = ((T(m^, s), o"(r?i^, s)) a diagonal matrix. 



(26) 



(27) 



In deriving (|2^) in [1^] one was concerned essentially with the right hand 
cut, responsible for unitarity in the corresponding channel. As a result, the 
imaginary part of the amplitudes to be used in (pO|), above the lightest 
threshold (in our case the vrvr one, s = 4m^) and below the highest one 
(s = Am\), was restricted to come only from unitarity, (p3|), neglecting the 
left hand cut contribution to the imaginary part that appears in Tcc,i and 
Tcn,i for s < — 4m^. 

However, we have maintained the left hand cut contribution to the imag- 
inary part of Tec, 4 and Tcn,A below the KK threshold. One way to see how 
large is the resulting deviation from unitarity is to check the value of the 
inelasticity in the energy region 4m^ < s < 4m|^ for (I,J)= (0,0) and (1,1) 
where two channels appear. In both cases the deviation from 1 is smaller 
than 1%. 

The rjT] intermediate state has been included in the 0{p'^) ChPT ampli- 
tudes. However, for the (0,0) channel for -^s > 2m,y this state gives further 
contribution to the imaginary part of our amplitudes in addition to the one 
expressed in (p3[). This means that (pO|) with only the tttt and KK states 



vrvr S-wave phase 



does not fulfill unitarity strictly for ^/s > 2m^ ~ 1.1 GeV. The infiuence 
of the rjr] state is particularly significative in the KK 
shifts, Fig. 6. We will come back to this point later. 
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Omiting the I, J labels, the relation between the T and the S matrix 
elements for a two channel process, in our case vrvr and KK for (I,J)=(0,0) 
and (1,1), is given by 



511= I + 2ia{ml,s)Tii 
5*22= I + 2ia{'mj^,s)T22 



Si2= 2i ^a{ml,s)a{ml,s)Ti2 (28) 
To accomplish unitarity the 2x2 S'-matrix can be written |2^ as 

From ( p8|) and ( p9|) we obtain the phase shifts for vrvr vrvr (^i) and vrvr 
KK {Si + Ja) for (I, J) = (0,0) and (1,1). For (I,J)=(2,0) only the vrvr 
channel is necessary, when omiting multipion states. In this case, it is enough 
to consider the first equality of (p8|) to obtain the phase shifts with 5ii = e*^*^. 



4. Fit, phase shifts and inelasticity 

In ChPT the experimental values for the Lj coefficients come from 0{p'^) 
fits to low energy experimental data. Here we fit the Lj constants to ex- 
periment in a much broader energy interval and with an expression valid to 
all orders. Hence, differences are expected between our fitted Li parameters 
and the values quoted from ChPT. 

Furthermore, our approach is not cross symmetric. This implies that 
contributions from the left hand cut of order higher than 0{p^) are effectively 
reabsorbed in the values of our Lj coefficients. This point has been studied 
in [^8| with the conclusion that the value of the Lj obtained from a non cross 
symmetric method are influenced by this reabsorption procedure of the left 
hand cut. In this way, the value we quote for our Lj constants has to be 
taken with care when comparing with the values of the Li from ChPT. 

We have used simultaneously the phase shifts of the vrvr vrvr with 1=0 
and 1, Figs. 4 and 5, to fit the value of our free parameters: Li, L2, -^3, L4, L5 
and 2Lq + Lg. The fit has been done using MINUIT. In the energy region 

= 500-950 MeV the data from different experiments for S-wave vrvr phase 
shifts are incompatible. Given that situation, we have taken as central value 
for each energy the mean value between the different experimental results 
[29-34]. For ^/s = 0.95-1 GeV, the mean value comes from 33|. In both 
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cases the error is the rriciximum between the experimental errors and the 
largest distance between the experimental points and the average value. 

The quoted errors in the value we have obtained for the Lj coefficients 
is just the statistical one. 

The fit is pretty good, as can be seen in Figs. 4 and 5, with = 1-3 
per degree of freedom. The values we obtain at the Mp scale and in units of 
10~^ are 



1 = U.^Z_Q02 
r _ 1 qc+0.02 

Lg = -3.24 ± 0.04 

L4 = 0.20 ± 0.10 

L5 = O.Olgi 

2Le + Ls = Omttil (30) 

The small errors for Li, L2 and L3 arc due to the strong constraints imposed 
by the small errors in the experimental data of the (5ii,7r7r phase shift, Fig. 4. 
The values of ChPT are 



Li = 


0.4 ± 0.3 


L2 = 


1.4 ± 0.3 


Lz = 


-3.5 ± 1.1 


L4 = 


-0.3 ± 0.5 


L5 = 


1.4 ± 0.5 


2Lq + Lg = 


0.5 ± 0.7 



(31) 

We can see that our values, taking errors into account, are compatible 
with those from ChPT, and then we can guarantee a good behaviour at low 
energies for our predictions. 

Using these values for Lj we also describe correctly phase shifts, scatter- 
ing lengths and form factors for the I = J = and I = J = 1 channels, as 
we will see later. 

In Fig. 4 wc show our fitted 5i for / = J = 1, from the two pion 
threshold up to 1.2 GeV and we see a good agreement with the experimental 
data which are dominated by the presence of the p(770) that we reproduce 
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nicely. In this channel we obtain that the influence of the KK channel 
coupled to vrvr is negligible. 




:iO0 400 000 600 700 800 900 1000 1100 1^00 



Ecn.[MeV] 

Figure 4: Phase shift for vrvr vrvr in / = J = 1. Data: 




300 400 000 600 700 600 000 1000 1100 1200 



Ecn>[MeV] 

Figure 5: Phase shift for vrvr — >■ vrvr in / = J = 0. Data: empty pentagon 



[p9| , empty circle |Q , full square , full triangle |Q , full circle represents 
the average explained above. 
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In Fig. 5 we show 6i for / = J = 0, also from two pion threshold to 1.2 
GeV. The agreement with experiment is quite satisfactory showing clearly 
the presence of the /o(980) resonance as a strong jump in the phase shift 
around 1 GeV. To get this resonance it is essential to include the kaons, and 
then to unitarize with coupled channel as we do. In this figure we also have 
plotted data from |2^, but since no error is quoted we have not included 
this data in the fit. 

Now, once we have fixed the Li from the fit we predict other magnitudes. 




1000 1020 1040 1060 1080 1100 1120 1140 1160 IIBO 1200 



E„„[MeV] 



Figure 6: Phase shift for vrvr — > KK in I = J = 0. Data: full square |35], 
full triangle p6| . 



In Fig. 6 the phase shift for the KK — > vrvr scattering, 6i + 62, is shown 
for I = J = 0. In this figure one sees clearly the r]r] threshold. This process 
is the most sensible to the rjT] intermediate state, contrary to what happens 
with the TTTT phase shifts. Fig. 5, where its inclusion is almost negligible. 
One way to realize the influence of this channel in the former phase shifts 
is to cancel the imaginary part to T4 coming from the intermediate rj?] s- 
channel loop for ^/s > 2m^ ~ 1.1 GeV. In this way the dashed line curve 
in Fig. 6 is obtained, which agrees very well with data. This is telling us 
that the inclusion of the rjr] channel in the unitarization procedure for the 
(0,0) channel is important for studying the KK — > vrvr scattering. In any 
case, as explained above, the threshold is very close to 1.2 GeV where other 
intermediate states, as four pions, are also important and should be included 
as well. Hence, we think that the inclusion of the rji] threshold in eq. (l2C 
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should be done when gomg to higher energy, that is, when extending the 
model for energies higher than 1.2 GeV. 

In Fig. 7, (1 — (??oo)^)/4 is shown. Our results display the same tendency 
as the experimental data, particularly when taking into account the large 
experimental errors. 




In Fig. 8 we show the vrvr phase shift with 1=2, J=0. The agreement with 
experimental data is fair. Contrary to the other channels we have shown, 
in this case no resonances appear and there is only the tttt channel. So, in 
this case our result (apart of differences on the Li values) is the same than 
in the lAM 0. 

We have also calculated the scattering lengths for the three channels 
unitarized in this work, (I,J)=(0,0), (1,1) and (2,0). We denote them by aj. 
In Table I we show the value we obtain for aj together with the experimental 
and the ChPT values to 0{p'^). We see in this table that a good agreement 
with experiment is accomplished. Our values are also close to the ones 
from ChPT as one should expect because for low energies we recover the 
chiral expansion. 
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Figure 8: Phase shift for vrvr vrvr in / = 2, J = 0. Data: cross |38|, empty 
square [39]. 



Table I: Comparison of scattering lengths in different channels 



«5 


ChPT 


Our results 


Experiment 


al 
«o 


0.20 ±0.01 
0.037 ± 0.002 
-0.041 ± 0.004 


0.210 ±0.002 
0.0356 ± 0.0008 
-0.040 ± 0.001 


0.26 ±0.05 
0.038 ± 0.002 
-0.028 ± 0.012 



5. Calculation of the scalar and vector pion form factors 

The scalar and vector form factors of the pion are defined respectively 



as 



and 



tt^{p)tt\p) out \m{uu + dd)\ O) = 6^^r{s) 



7r*(p')7r'(p) out 



97, y U 







(32) 



(33) 



with m = {niu ± md)/2 and e^^^ the total antisymmetric tensor with three 
indices. 

Assuming elastic unitarity (valid up to the KK threshold and neglecting 



multipion states) and making use of the Watson final state theorem |21| the 
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phase of r(s) and Fy^s) is fixed to be the one of the corresponding partial 
wave strong amphtude: 



Imr(s + ie) 
Im Fv{s + ie) 

The sohition of (|3^) is well known 




= tan(5[;Rer(s) 

= tanjj ReFy(s) (34) 
and corresponds to the Omnes type 



With 



T{s)= Po{s)nois) 

Fv{s)= Pi{s)ni{s) (35) 



aO^=exp -/ —^J^-±-\ 36 

vr Uml s' s' - s - te\ 



In ( |35| ) ^'o('S) and -Pi(s) are polynomials of degree fixed by the number of 
subtractions done in ln{f]o('S)} and ln{$7i(s)} minus one, and the zeros of 
Fv and T. For n = 1, Pi{s) = 1. This follows from the normalization 
requirement that r(0) = -Fy(O) = 1 and the absence of zeros for those 
quantities. 

Since we have a prediction for the phase shifts we can calculate the 
dispersion integral (|3^ ) and obtain the pion form factors for both the scalar 
and vector cases. The results are shown in Figs. 9 and 10. The Omnes 
solution assumes the phase of the form factor to be that of the scattering 
amplitude, and that is true exactly only until the first inelastic threshold. 
The first inelastic threshold is the 47r one. However, as it was already said, 
its influence, in a first approach, is negligible. The first important inelastic 
threshold is the KK one around 1 GeV. This is essential in / = J = but 
negligible in / = J = 1. This inelastic threshold, as discussed above, has 
been included in our approach and it is responsible for the appearance of 
the /o(980) resonance, as it is clearly seen in Figs. 5 and 10. 

Due to the presence of this important KK threshold, overall for the 
I = J = sector, eqs. (^5|) are strictly correct up to -/s = 2mx, which 
is indicated as a dashed-dotted vertical line in Figs. 9 and 10. Up to this 
energy, we see the clear appearance in Figs. 9 and 10 of the /o(770) and 
/o(980) respectively. In the case of the Fv{s) the agreement with existing 
data is quite satisfactory. Above the KK threshold one expects deviations 



from (35) due to the opening of this inelastic channel. However, for the 



17 



vector form factor we still see a rather good agreement with data and the 
deviation should be adscribed to the presence of the (J resonance above 
1.2 GeV. On the other hand, the result obtained for the vector form factor 



is similar to the one recently obtained in |17| using another phase shift 
expression, taking into account possible uncertainties coming from orders 
higher than "p^ in ChPT. For the I = J = channel the most dramatic 
influence of the opening of the KK threshold is the appearance of the /o(980) 
resonance, what happens a little below KK threshold. Hence its appearance 
in Fig. 10 is well accommodated in our assumption of elastic unitarity for 
the Watson theorem, which we have used to evaluate the form factors. So 
that we do not expect large deviations from our results even above the KK 
threshold up to the appearance of new / = J = 0, /o, resonances higher in 
energy, typically around ^/s ~ 1.3 GeV. In Fig. 10 the dashed line represents 
the scalar form factor unitarizing only with pions to obtain the (Joo.tttt phase 
shift, in the line of the works |l^, 27 1 and we see a very large influence of 
the KK channel through the /o resonance which is even substantial around 
\^/s\ = 500 MeV. 




-400 -200 200 400 600 BOO 1000 l^dOO 1400 1600 



Figure 9: Vector pion form factor. The vertical line shows the opening of 



the i^/C threshold. Data: |0|. 
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Figure 10: Scalar form factor. The dashed curve is the result unitarizing 
only with pions. The solid line is the full result with both pious and kaons 
in the intermediate state. The vertical line shows the opening of the KK 
threshold. 



6. Conclusions 

In this work we have calculated the KK scattering amplitude to next to 
leading order in ChPT. We have seen that from the KK threshold, due to 
the large kaon mass, very large corrections appear. This in principle should 
imply that this O(p^) calculation is unlikely to be useful. However, we see 
that this is not true and that one can obtain accurate results when the ChPT 
calculations are supplied with some suitable non perturbative unitarization 
scheme as the one used here. In this way, we have successfully described 
the TTTT — TTTT phasc shifts for (I,J)=(0,0), (1,1) and also the phase shift 
for TTTT KK in (I,J) = (0,0) and the inelasticity in good agreement with 
experiment up to ^/s = 1.2 GeV. The scattering lengths for tttt with 1=0,1 
and 2 are also calculated in agreement with experiment and former ChPT 
calculations. We have finally computed the scalar and vector form factor of 
the pion making use of the above calculated vrvr phase shifts in an Omnes 
representation and the results have also been satisfactory. The scalar form 
factor strongly shows that the KK channel is essential in order to reproduce 
the /o resonance. 
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